When can preheating affect the CMB? 
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We discuss the principles governing the selection of inflationary models for which preheating can 
affect the CMB. This is a (fairly small) subset of those models which have non-negligible en- 
tropy/isocurvature perturbations on large scales during inflation. We study new models which belong 
to this class - two-field infiation with negative nonminimal coupling and hybrid/double/supernatural 
inflation models where the tachyonic growth of entropy perturbations can lead to the variation of 
the curvature perturbation, TZ, on super-Hubble scales. Finally we present evidence against recent 
claims for the variation of TZ in the absence of substantial super-Hubble entropy perturbations. 



I. INTRODUCTION 

The cosmological landscape is now dominated by a 
myriad of inflationary models [0, each with slightly dif- 
ferent genetics but a common origin, forged in the heyday 
of Grand Unifled field theories ||]. Inflationary models 
are flexible, hardy and spawn & multiply with remarkable 
facility - features flowing from their scalar field DNA. 

And while the recent M-theory and brane-world rev- 
olution has again stimulated interest in alternatives to 
inflation |^ it remains true that the simplest infiation- 
ary models provide a very good fit to current cosmo- 
logical observations Q]. Indeed, the Cosmic Microwave 
Background (CMB) and large scale structure observa- 
tions presently show no particular signature - such as pri- 
mordial non-Gaussianity - which might allow us to single 
out a particular model or falsify inflation as a paradigm. 
It is therefore good science to seek new ways to select the 
"fittest" of inflationary models. 

Preheating after inflation may provide exactly such 
a selection rule. Preheating is probably the most vio- 
lent of putative phases in cosmic history It re- 
quires non-perturbative quantum field theory pushed to 
its non-equilibrium limits 0. The exponentially rapid 
population of sub- and super-Hubble wavelength modes 
makes it possible to produce very massive particles in 
large quantities and hence over-produce dangerous relics 
such as modulini or gravitini In addition there is 
some (admittedly speculative) evidence that the induced 
parametric growth of small-wavelength metric perturba- 
tions around the Hubble scale may lead generically to 
runaway production of primordial black holes (PHB) 
If correct this will rule out wide regions of the preheating 
parameter space. 

At the other end of the cosmic scale, a natural question 
is whether preheating can affect the modes of the metric 
perturbations which induce the intricate anisotropies in 
the CMB, and if so, what are the criteria? In this paper 
we will discuss the conditions where preheating can lead 
to the growth of metric perturbations on cosmic scales. 

The answer to this question appears to be intimately 



linked to the issue of correlations between adiabatic and 
entropy perturbations from multi-field inflation. 



II. GENERAL PRINCIPLES 

In this section we present the conditions known to be 
sufficient for preheating to affect the CMB; distilled from 
the recent literature jlOj- |^l[. Perhaps the best way to 
discuss the impact on the CMB is via quantities which 
are time-independent in the presence of only adiabatic 
perturbations on large scales. Such quantities are used 
to normalize the inflationary models to the large-angle 
COBE measurements of the CMB. 

Two common choices for such quantities are Q, the cur- 
vature perturbation in the constant density hypersurfaces 
{5p = 0), and the comoving curvature perturbation TZ. 
These two coincide in the fc — ^ limit ||2^ (modulo a 
minus sign) and we will use TZ. 

The existence of entropy perturbations is crucial for 
preheating to affect the CMB, since in general p3] 
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(2.1) 



where the above limit is understood as fc — > 0; H ^ p, p are 
the Hubble rate, the total pressure and energy density 
respectively and S is the total (dimensionless) entropy 
perturbation deflned by 



5 = i/('^-^ 
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Now in the case of two minimally coupled scalar fields 
tfi and tf2^ this can be given explicitly ||22| as 
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where a and s are "adiabatic" and "entropy" fields, de- 
fined by 
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da = {cos0)dLpi + {sm0)dLp2, 
ds — ~'{sm9)d(pi + {cosd)d(p2 



(2.4) 
(2.5) 



Here 9 is the angle of the trajectory in (0, x) held space, 
satisfying tan 9 — (^2 /^i • Eq- ( ^.3| ) means that variation 
of TZ requires not only a large scale entropy perturbation 
ds, but also a non-straight trajectory in field space. In 
the minimally coupled two-field system with an effective 
potential V{ipi,(p2), the Fourier modes for the adiabatic 
and entropy field perturbations satisfy [E2| 
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where k is a comoving momentum, a is a scale factor, G 
is a Newton's gravitational constant, and 

K<T = (cos^ 9)V^,^, + (sin 29)V^,^, + (sin^ 9)V^,^, , (2.8) 
Vss = (sin^ e)V^,^, - (sin 29)V^,^, + (cos^ 0)V^,^, . (2.9) 

$ is a gravitational potential in the longitudinal gauge, 
satisfying 



$ + iJ$ = ATTGa6a. 



(2.10) 



Eq. ( 2.1C ) indicates that the gravitational potential is 
sourced by adiabatic field perturbations. 

When the effective mass of 6s is light relative to the 
Hubble rate H — d/a during inflation, i.e.. 
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the entropy field perturbation is not suppressed on super- 
Hubble scales during inflation. Then during preheating 
if 6s is resonantly amplified due to a time-dependent ef- 
fective mass, this can lead to the growth of TZ on large 
scales, thereby altering the power spectrum normaliza- 
tion or even leaving the model incompatible with the 
large-angle CMB. Note that the adiabatic field perturba- 
tion is sourced by the entropy field perturbation, thereby 
stimulating the growth of $ through Eq. ( 2.10 ). In con- 
trast, if the entropy perturbation is heavy during infla- 
tion, (^5 > ij2) then \6s\ ^ a"^/^ and the growth during 
preheating means that the change of TZ is negligible be- 
fore backreaction ends the resonance. 

In Sec. IV we present new classes of models which 



have strong preheating or tachyonic growth but simul- 
taneously have a light entropy perturbation in the pre- 
ceding inflationary phase. To begin with we shall analyze 
the evolution of TZ in the massive chaotic inflationary sce- 
nario in the next section. 



III. AN EVALUATION OF CLAIMS FOR 
VARYING 7^ IN THE ABSENCE OF ENTROPY 
PERTURBATIONS 

There have been recent claims by Henriques and Moor- 
house (HM) |]l^ that TZ or ^ will vary during reheating or 
preheating even in the absence of large-scale entropy per- 
turbations when going beyond linear perturbation theory 
and taking into account the quantum-to-classical transi- 
tion. 

If correct this would have a profound impact on infla- 
tionary cosmology Our aim in this section is to 
evaluate these claims. To do so let us consider metric 
preheating in the massive chaotic inflationary scenario 
with a standard four-leg interaction: 
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Strong amplification of the x fluctuation requires a reso- 
nance parameter q = (fp' / {Am^) 1 at the beginning of 
preheating . In this case the field x is heavy during in- 
flation relative to the Hubble rate , which means that 
X is strongly suppressed during inflation {x ^ a^^/"^), 
thereby leading to ^ ~ and {Vss) — {Vxx) — 9^4''^ ^ 
H^. Therefore large scale entropy perturbations are ex- 
ponentially suppressed during inflation, which safeguards 
nonadiabatic growth of super-Hubble curvature pertur- 
bations, as found by Eq. (^.3|). Hence, at the linear level, 
7?. ~ to high precision in this model. 
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FIG. 1. The evolution of a super-Hubble curvature pertur- 
bation iZ = k'^^^TZ and the gravitational potential k^^^'^ 
during inflation and preheating in the model (3.1) with 
g = 5.0 X 10~ . The initial conditions are chosen to be 
(j) = 3Mpi and x = 10~"^Mpi, in which case inflation ends 
around mt ~ 20. Inset: The plot of TZ for large initial vari- 
ance, {54)'^) = 10~*(^o- 



Nevertheless we have to caution that TZ in Eq. (2.3) is 
the result of the first order perturbation theory. If we in- 
clude second order backreaction effects in the background 
evolution equations, these give rise to second order terms 
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such as (50^) in the time derivative of TZ Since these 
variances are not suppressed during inflation (due to the 
short-wavelength contributions), they may provide addi- 
tional source terms for the curvature perturbation. This, 
in essence, is the origin of the claims by HM. 

To check these claims we performed our numerical 
simulations implementing the second order backreaction 
effects as spatial averages in the evolution equations 
P,[7| p9|j20[ | and compared the evolution of the cosmolog- 
ical perturbations found using different (but equivalent 
at linear order) equations of motion. 

In Fig. 1^ we plot the evolution of a super-Hubble 
curvature perturbation using its definition, i.e., TZ — 
$ — H/H{^ + H^). The variance (Sx^) grows by para- 
metric resonance until backreaction becomes important, 
while the inflaton fluctuation is not excited unless rescat- 
tering is taken into account [^ . 

The coherent oscillation of the inflaton condensate is 
destroyed once backreaction starts to dominate. This can 
affect the evolution of H, $ through Eq. (2.10) and also 
TZ. In fact the curvature perturbations exhibits small 
oscillations, but it remains roughly co nsta nt (see Fig. |^). 
When we use the first order equation ( ^.3| ), we found that 
TZ is conserved without oscillations. 

This is in stark contrast to the simulations of HM jT^ ] 
who found that TZ varies after the energy density of the 
inflaton condensate drops below its variance {ScjP). Ac- 
cording to their numerical results using similar second 
order approximations such as ours, the decrease of TZ oc- 
curs in the preheating stage. 

A large difference between the two investigations is 
that the initial energy density of the variance {5<lP) is 
smaller than <jy^ only by four orders of magnitude in 
Ref. while we argue that the typical size is regu- 

lated to be (50^) ~ ~ IQ-i^M^j |24), which is by ten 
orders of magnitude smaller than 0q . 

Since (50^) does not exceed 0q during preheating in our 
simulations, we do not flnd the decrease of TZ claimed by 
prf . Using the initial conditions of HM we checked that 
curvature perturbations exhibit small decrease after 
drops under {Scj)^) (see the inset of Fig. 1), but we still 
do not flnd the extensive changes reported in [ pT[ . 

We did flnd that implementing the numerics is subtle 
and subject to artiflcial instabilities. Further, not all 
deflnitions or equations for TZ are equally suitable for 
numerical implementation, some being more susceptible 
to these instabilities. We therefore argue that TZ probably 
does not evolve on large scales in the absence of entropy 
perturbations even at second order. This is consistent 
with the standard view based on causality 0. 



*The standard view simply associates t he sec ond order terms 
to radiation with a temperature T oc (Sx^)- In the absence 
of large scale entropy perturbations in this radiation fluid no 
changes are induced in TZ. 



We end with one caveat, however: full lattice simula- 
tions show that the growth of 5x leads to the excitation of 
inflaton fluctuations via rescattering, thereby satisfying 
the condition {5(j)^) > 0q around the end of preheating 
. It would be worth investigating further whether any 
change in TZ occurs in such a situation, although if they 
do, they are likely to be small. 



IV. MODELS WITH THE CMB AFFECTED BY 
PREHEATING 

A. Chaotic inflation with self-interaction and 
nonminimal coupling 

Achieving a light entropy perturbation during infla- 
tion is typically rather difflcult in chaotic inflation mod- 
els. However, this picture is modiflcd if one takes into 
account nonminimal coupling [ p^ . Let us consider the 
quartic chaotic inflationary scenario in the presence of a 
nonminimally coupled scalar field x coupled to 4> with 
effective potential: 



In this case the effective mass of x is given by 
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where we used the approximation R w 12H^ which as- 
sumes X ~ s-iid is valid to zero order in the slow-roll 
parameters. Since the term decreases faster than the 
g'^cfP term, it is possible for x to be light relative to H 
during inflation by allowing negative values of ^. When 
^ = it was shown in Rcfs. |l6yi§| -|20|| that super-Hubble 
cosmological perturbations probed by CMB experiments 
can be amplifled around the center of the flrst resonance 
band, g'^/X = 2 (see also Ref. (25)). For g^/X > 8, us- 
ing the Hartree approximation, the growth of sub-Hubble 
field perturbations shuts off the resonance before super- 
Hubble metric perturbations are enhanced |po| , p6| . 

However, this picture is modified by a negative non- 
minimal coupling for x, which makes it possible to avoid 
the inflationary suppression of the entropy perturbation 
even for g'^/X > 8. For example, when g^/X — 18 and 
^ = —0.12 shown in Fig. ^, the amplitude of the supcr- 
Hubblc Sxk mode at the end of inflation is larger than in 
the C = case by about 10 orders of magnitude. In this 
case, large-scale curvature perturbations exhibit nonadi- 
abatic growth after Sx catches up (50 (see Fig. ||). 

When g^ / X 3> 1, rather large negative nonminimal 
coupling (^ < — 1) is required to make the x mass light. 
In this case, if the ^R term gives almost the same contri- 
bution as g^(j>^ at the beginning of inflation, it is impossi- 
ble to avoid the suppression of x due to the decrease of £,R 
compared to g^(j)^. When g^0^ < —^R initially, strong 
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amplif ication of TZ occurs from the beginning of inflation 
1^,^. Therefore when g^/X ^ 1 with large negative 
nonminimal coupling, the growth of TZ during inflation 
is significant rather than during preheating. The density 
plot of the parameter range where curvature perturba- 
tions grow during inflation and preheating is shown in 
Fig. I 
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FIG. 2. The evolution of super-Hubble cosmological per- 
turbations for g^/X = 18 and ^ = —0.12 with initial conditions 
(j> = 4Mpi and x = 10~^Mpi. Nonminimal coupling makes the 
X mass light, thereby leading to the amplification of 7?. during 
preheating. Here x is the dimensionless conformal time. 



In the case of the massive inflaton and a g^c^^x^ cou- 
pling, the effective mass of x is approximately given by 
miff « [g^ -I- 167r^(m/Mp;)^](^^ during inflation. This im- 
plies that it is not possible to achieve both the require- 
ment of a light X niass during inflation and a large x niass 
(g 3> 1) after inflation and hence TZ is not enhanced dur- 
ing preheating even in the presence of nonminimal cou- 
pling. 

This situation changes, however, if we consider a cubic 
coupling of the form g'^4'x'^ instead. This still leads to 
preheating but now we can make rn^g light during infla- 
tion in certain regions of the {g/m^S^) parameter space, 
analagous to those in Fig. |^. 



B. Hybrid and Double inflation 

In hybrid and supernatural inflationary models 
|p9|-^, a symmetry breaking transition occurs in the 
presence of a second scalar field, N . This leads to a neg- 
ative coupling/tachyonic instability in which the longest 
wavelengths grow the strongest. Hence, even in the case 
when the spectrum of N fluctuations is blue, this long- 
wavelength instability may lead to variation of TZ in cer- 
tain cases. 

To investigate this in detail, consider Linde's hybrid 
inflation potential |29| with two scalar fields and N: 
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FIG. 3. A density plot in the {g^ /X, 5)-plane where the su- 
per-Hubble curvature perturbations TZ is enhanced or not. 
The black color and the dark gray colours correspond to 
ranges where TZ grows during inflation, and during preheating, 
respectively. The region where TZ does not vary is shown by 
the light gray color. Note that there exist parameter ranges 
around the center of resonance bands {g^ /X — 2n^) where TZ 
grows during preheating even for g^ /X > 8. However, these 
regions do not persist for g^/X > 40 as explained in the text. 



^'"^) +^5^^''^' + ^-^. (4.3) 



Note that the supersymmetric version of the hybrid in- 
flation corresponds to the case, g^/A — 2 |30j. Hereafter 
we shall analyze the case where g^ and A are the same 
order. 

Inflation takes place due to the slow-roll evolution 
of (j) before the field reaches a bifurcation point 0c = 
M/g. If the condition m^^^ M^jX is satisfied, the 
H ubble co nstant at = is given by i7 = iJo = 
V27r/(3A)MVMpi. It is convenient to normalize the 
masses of two fields </) and relative to as 



m2 _ iXrn^Ml, 



/3: 



M2 



3AM| 



(4.4) 



a is required to be smaller than unity in order to lead to 
sufficient inflation for (j) > (j>c- The evolution of the fluc- 
tuation 5N depends on the value of (3 as we will classify 
below. 



1. Light N field (Double inflation) 

When the N mass is light (/3 < 1), the field TV and its 
large scale perturbation are free from inflationary sup- 
pression {(j) > 4>c)- In this case the coupling A is con- 
strained to be small, A < (M/Mpi)^. In order to end 
inflation due to a rapid rolling of the field TV after the 
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symmetry breaking, one requires the "water-fall" condi- 
tion < AmM^j !§. If we adopt the COBE normal- 
ization, gM^/{X^/^m'^M^{) ~ 3.5 x 10"^ which comes 
from the single-field fluctuation (50, the water-fall condi- 
tion is not typically satisfied due to small A. In this case 
the second stage of infiation occurs after the symmetry 
breaking (double inflation). 
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the oscillating phase. This illustrates the importance of 
the spinodal instability region, during which TZ and are 
amplified as well as field perturbations. 



2. Heavy N field (Hybrid inflation) 

When /3 ^ 1, the field N is exponentially suppressed 
during inflation > (f>c)- In this case the water-fall 
condition is typically satisfied after symmetry breaking, 
which corresponds to the original hybrid inflationary sce- 
nario p9| where inflation ends due to the rapid rolling of 
the fleld N. 

The evolution of fluctuations for (p < (pc can be de- 
scribed similarly as the models of spontaneous symmetry 
breaking analyzed in Ref. [Q. The field N has practi- 
cally no homogeneous component a,t (p = (pc- In this case 
the decomposition of N between the homogenous field 
N{t) and the fluctuational part 6N{x, t) is not necessar- 
ily valid. Rather one is required to go beyond perturba- 
tion theory by considering the full spatial distribution of 
the field ^(x,t). 



FIG. 4. The evolution of super-Hubble cosmological per- 
turbations in the case of a light A'' mass with potential ( [l.3| ) 
as a function of i = \/2Mt using the Hartree approxima- 
tion. We choose the initial conditions to he (p = 1.5(/f>c and 
A'^ = M / \/\. We also show the evolution of adiabatic 
and entropy field perturbations, denoted by 5a and 5s, re- 
spectively. 

Let us present a concrete example. In Fig. ^ we plot 
the evolution of Tl and $ for M = 5.0 x 10"'^Mpi, g"^ = 
A — 8.0 X 10~^^, and m — 0.2M, corresponding to mass 
parameters a = 0.07 and /3 = 1.5. In this case the first 
stage of inflation {(p > (pc) continues around i = \f2Mt ~ 
33, after which the second phase of inflation takes place 
until i ~ 130. The total number of e- folds is ~ 70 for 
the initial value, (p> = 1.50c- 

In Fig. ^ we find that the entropy perturbation Ss is 
not strongly suppressed and begins to increase after the 
symmetry breaking due to tachyonic growth of the fluc- 
tuation SN. This can lead to the amplification of super- 
Hubble curvature perturbations until the backreaction of 
field fluctuations shuts off its growth. 

Note that we implement the field backreaction effect 
at second order in our simulations. Since the fluctuation 
SN is efficiently amplified in the tachyonic instability re- 
gion, its growth terminates before the field TV reaches the 
potential minimum. The excitation of entropy pertur- 
bations also leads to the amplification of adiabatic field 
perturbations da since they are correlated each other. In 
Fig. ^ we clearly find that the gravitational potential ex- 
hibits nonadiabatic growth sourced by the growth of Sa. 

During reheating no enhancement of super-Hubble cos- 
mological perturbations occurs in this case, since field 
perturbations are already excited until the beginning of 
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FIG. 5. The evolution of the large scale curvature per- 
turbation in the heavy A'' case with M = 3.0 X 10"®Mpi, 

= X — 1.0 X 10~* using the Hartree approximation to eval- 
uate the bac kreaction. The result (i) is obtained by integrat- 
ing Eq. (2.3), which is different from the result (ii) obtained 
by the definition 7^ = $ - H/H( ^ + H $). In the case (ii) the 
field a is strongly dumped in Eq. (2.10) due to overestimation 
of the field variance, thereby leading to the decrease of "I> and 
TZ. In the case (i) the damping of a provides a source term 
for TZ. 



If the perturbative approach based on the Hartree ap- 
proximation is applied, the effective mass of the field A^ 
around = is approximately given by mj^ ~ 3A((5A^^) — 
Af^. Naively one may expect that symmetry can be re- 
stored due to the growth of SN with {SN"^) ~ M^/X, 
i.e., ~ 2M^. Numerically we have found that the 
Hartree approximation leads to symmetry restoration. 
However we have to caution that this approach neglects 
the rescattering effect and the formation of topological 
defects. The authors in Ref. fS^ found that the field dis- 
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tribution reaches the potential minimum at iV = M/TA 
and (j) = instead of restoring symmetry when using the 
3-D lattice simulations without metric perturbations. It 
was also pointed out that tachyonic growth of the field 
fluctuation is typically followed by only a single oscilla- 
tion of the field. 

In the Hartree approximation, the backreaction due 
to the growth of long wavelength fluctuations prevents 
the field from reaching the potential minimum due to an 
overestimation of the field variance. In Fig. ^ we find that 
the evolu tion of TZ differs significantly when we integrate 
Eq. (2.3) or evaluate TZ in terms of $. This disagreement 
obviously shows the limitations of linearized gravity using 
the Hartree approximation. It is certainly of interest to 
analyze the evolution of super-Hubble cosmological per- 
turbations in the heavy N case including full nonlinear 
effects along the lines in Ref. . 



V. CONCLUSION AND DISCUSSION 



tion such as discussed in Sec. IV B - a natural ques- 
tion arises: "what are the correct perturbed Einstein 
field equations?" It would be very useful to develop a 
formalism in which these quantum fields were treated 
non-perturbatively and self-consistently. 

Finally, does inclusion of rescattering leads to extra 
variation of TZ or does causality protects TZ (and suitable 
nonlinear generalizations) to all orders in perturbation 
theory |35j in the absence of large-scale entropy pertur- 
bations? These interesting issues are left to future work. 
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We discussed general principles governing the condi- 
tions under which the large scale curvature perturbation, 
TZ - used to normalize models to the large angle CMB 
anisotropics - is amplified after (and during) inflation. 

This is closely allied with the evolution of en- 
tropy/isocurvature field perturbations Ss during the in- 
flationary phase. If the mass of 5s is light relative to 
the Hubble rate, it sources the variation of the curvature 
perturbation when it is amplified during preheating or in 
the tachyonic instability region. 

In the model with V = ^rn?(jP + ^ff^^^X^j super- 
Hubble entropy {5s) modes are exponentially suppressed 
during inflation for the coupling g required to have ef- 
ficient preheating. Therefore growth of 5s during pre- 
heating is not sufficient to lead to the variation of TZ on 
large scales. This situation is not altered in the model 
V = \w?(lP' + \g'^(t>^X^ + g'^^'^X where the fields evolve 
along the straight line {9 ~ 0) during inflation |2|] . 

In the model with V = iA0* + \g'^(tP'X^ + i^^X^ we 
found that parameter ranges where large scale curvature 
perturbations are amplified during preheating are signif- 
icantly wider for negative nonminimal couplings. 

We also studied hybrid/double infiationary models 
where tachyonic/spinodal instability of the symmetry 
breaking field N leads to the growth of entropy field per- 
turbations. When the field N is light in the first stage of 
inflation, super-Hubble curvature perturbations typically 
exhibit tachyonic growth after the symmetry breaking. If 
the field N is heavy, which is the original version of the 
hybrid inflationary scenario, we flnd that linearized grav- 
ity theory using the Hartree approximation shows some 
limitation in the sense that the evolution of TZ is different 
depending on which deflnition of TZ is used. 

In the case of a purely quantum field with no classi- 
cal vacuum expectation value (vev) - such as a fermion 
field or a scalar field whose vev vanishes during infla- 
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